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ABSTRACT 

We examine the relation between the density variance and the mean-square Mach number in supersonic, 
isothermal turbulence, assumed in several recent analytic models of the star formation process. From a series 
of calculations of supersonic, hydrodynamic turbulence driven using purely solenoidal Fourier modes, we find 
that the 'standard' relationship between the variance in the log of density and the Mach number squared, i.e., 
0\ np /p — In (l + b 2 A4 2 ), with b = 1/3 is a good fit to the numerical results in the supersonic regime up to 
at least Mach 20, similar to previous determinations at lower Mach numbers. While direct measurements of 
the variance in linear density are found to be severely underestimated by finite resolution effects, it is possible 
to infer the linear density variance via the assumption of log-normality in the Probability Distribution Func- 
tion. The inferred relationship with Mach number, consistent with a p /p k, bM. with b = 1/3, is, however, 
significantly shallower than observational determinations of the relationship in the Taurus Molecular Cloud 
and IC5146 (both consistent with b » 0.5), implying that additional physics such as gravity is important in 
these clouds and/or that turbulent driving in the ISM contains a significant compressive component. Magnetic 
fields are not found to change this picture significantly, in general reducing the measured variances and thus 
worsening the discrepancy with observations. 

Subject headings: turbulence — ISM: structure — hydrodynamics — stars: formation — magnetohydrody- 
namics (MHD) — shock waves 



1. introduction 

The last few years have seen an increasing number of 
analytic models of the star formation process that use the 
log-normal density probability distribution function (PDF) 
produced by supersonic turbulent flows to predict statisti- 
cal quantities such as the initial and/or core mass function 
(e.g.lPadoan & Nordlundl2002tlHennebelle & Chabrierl2008 , 
2009) and the star formati on rate dKrumh olz & McKeej |2005 : 
Padoan & Nordlund 2009). A key assumption in these models 
is a relationship identified in early numerical studies between 
the PDF width - the density variance or standard deviation - 
and the Root Mean Square (RMS) Mach number A4 in su- 
personic, isothermal turbulence. The relationship is generally 
assumed to be linear in the standard deviation of linear den- 
sity, i.e., 

a p /p=bM, (1) 

where b is a constant of order unity and density is scaled in 
terms of the mean, p. For a log-normal distribution, this is 
equivalent to 

a 2 = In (1 + b 2 M 2 ) , (2) 

where s = \n(p/p), such that a s is the standard deviation in 
the logarithm of density. 
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Apart from the ear ly empirical findings of 
IVazauez-Semadenil d 19941). iPadoan. Nordlund. & Jones! 
d 1997b!) and iPasso t & Vazquez-Sem adenil d 19981) . there 
is no clear reason why the relationship should be of this 
form. Mathematically, the appearance of a log-normal 
distribution can be understood as a consequence of the 
multiplicative central limit theorem assuming that indi- 
vidual density pertu r bation s a re independent and random 
( Vazquez-Sema denil 11994]: | Passot & Vazque z~-SemadenH 
119981: iNordlund & Padoanl 1 19991) ! In physical terms this 
has been interpreted as meaning that density fluctuations 
at a given location are constructed by successive passages 
of shocks wit h a jump amplitude independent of the local 
density (e.g. [B allesteros-Par edes et alj 120071 iRritsuk et al.l 
120071: [Federrath et al.ll2.010t) . However it has not so far proved 
possible to analyti cally predict the relationshi p based on these 
ideas (though see Pado an & Nordlund! 120091) . Thus, a com- 
mon approach in numerical studies of turbulence — usually 
at a fixed Mach number — has been to measure the parameter 
b, assuming Eq. (|2), that gives best fitting log-normal to the 
time averaged PDF. However, reported estimates fo r b are 
widely discrepant. For example, Padoan et al. ( 1997b) found 
b ~ 0.5 while more recently Kritsu ket al.1 (!2007[) (at Mach 6) 
find a much lower valu e of b ~ 0.26 andlBeetz et al.l (12008 ) 
find b ~ 0.37, while IPassot & Vazquez-Semadenil (1998) 
f ound b ~ 1 (though with some confusio n ov er a s vs. Vp/p)- 

IFederrath. Klessen. & Schmidt! ([2008 ) and IFederrath et al.1 
(l2010r reconcile these results in part by the finding that the 
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width of the PDF depends not only on the RMS Mach number 
but also on the relative degree of compressible and solenoidal 
modes in the turbulence forcing, with 6=1/3 appropri- 
ate for purely solenoidal and 6 = 1 for purely compres- 
sive forcing. This is in keeping w ith e arlier discussions by 
Passot & Vazquez-Semadeni ( 1998) and Nordlund & Padoan 
d 19991) . the latter authors noting that "for compressional forc- 
ing at low Mach numbers (leading to an ensemble of sound 
waves), the standard deviation is expected to be equal to the 
RMS Mach number itself". 

Observationally, the log-normality of the 3D PDF is re- 
flected in the 2D column density PD F, for example as mea- 
sured from dust extinction maps (e .g. lLombardi et al J 120061 
120081 120 lOtlKainulainen et al.l2009l) - at least in earlier stages 
of molecular cloud evolution, suggesting that this phase could 
be dominated by roughly isothermal turbulence in which 
self-gravity is relatively unimportant. Onl y for seemingly 
more e volved clouds (including Taurus) do Kainulaine n et al.1 
(2009) see signifi cant tails at higher ( column) densities (sim- 
ilarly found by Lombard i et alJl20Tu) . However, measure- 
ments of the projected 2D variance (or PDF) cannot be di- 
rectly used to constrain the relationship with M ach num- 
ber. Recently, Bru nt. Federrath. & Pric e (2010a) (hereafter 
BFP) have shown how projection effects can be overcome 
to infer the 3D density variance from column density ob- 
servations, in turn leading to a method for extracting the 
unprotected (3D) density PDF fr om the observational data 
(iBrunt. Federrath. &" Price 2010b). This enables the relation- 
ship between the standard deviation in linear density and 
Mach number to be tested observationally, with initial ap- 
plication to Taurus finding b = 0.48 + ° dBruntl 120101). A 
similar method was employed by Padoan, Jones, & Nordlund 
( 1997a) to infer the 3D density variance from extinction mea- 
surements in IC5146, similarly finding 6 w 0.5. 

The problem with all of the above is that calculations — 
or observations — performed at a single (RMS) Mach num- 
ber can only ever assume the relationship given by Eqs. (O or 
(fl]l and cannot be used to co nstrain it unless a range o f Mach 
numbers are studied. Indeed, Lem aster & Stond (12008!) (here- 
after|LS08) — performing a series of calculations with Mach 
numbers in the range 1.2 < M < 6.8 — find a relationship 



a 2 = 0.72 In (1 + Q.bM 2 ) - 0.20, 



(3) 



based on a fit to measurements of the mean in the logarithm 
of density, s, as a function of A4, which we have here con- 
verted to a a s -M. relation using s = —a 2 /2. However, this 
three-parameter fit is clearly not unique, and it remains to be 
determined whether this or a similar relationship continues to 
hold at higher Mach numbers. 

The present study is motivated by a need to compare the 
theoretical predictions with the observational constraints. In 
particular LS08 only perform calculations up to M. w 6.8 
- corresponding to ID line Full-Width-Half-Maximum of 
~ 1.6 km/s (at 10K), which is rather low in terms of what 
is found in the real interstellar medium. In particular Tau- 
rus has M ~ 17, so a study going up to (at least) Mach 20 
or so is needed. Our aim in this paper is precisely this: To 
pin down the theoretical relationship - with as few assump- 
tions as possible - up to sufficiently high Mach numbers that a 
meaningful comparison can be made with observed molecular 
clou ds. Whilst additi onal physics such as non-isothermality 
(e.g. lScalo et alJ]T998l) . the multi phase nature of the interstel- 
lar medium and self-gravity (e.g. Klessen 2000; K ritsuk et al.l 



2010) are all expected to change the theoretical predictions 
at some level, the isothermal, non-self-gravitating case is an 
important reference point that remains theoretically uncer- 
tain. Furthermore, a clear prediction for this simple case can 
be used to gauge the relative importance of such additional 
physics in observed clouds. 

2. METHODS 
2.1. Log-normal distributions 
The log-normal distribution is given by 



p(s)ds 



1 



■. exp 



ds, 



(4) 



where s = \n(p/p) such that s and a s denote the mean and 
standard deviation in the logarithm of (scaled) density, respec- 
tively, and p is the mean in the linear density. The mean and 
variances in a log-normal distribution are related by 



and 



'p/p 



exp I a s 



(5) 
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2.2. Numerical simulations 

We have performed a series of calculations of supersonic 
turbulence, solving the equations of compressible hydrody- 
namics using an isothermal equation of state (with sound 
speed c s — 1) and periodic boundary conditions in the three- 
dimensional domain x,y,z e [0, 1]. Initial conditions were 
a uniform density medium p = p = 1 with zero initial ve- 
locities. Turbulence was produced by adding a random, cor- 
related stirring force, driving the few largest Fourier modes 
1 < k < 3 with a random forcing pattern, slowly changed 
according to an Ornstein-Uhlenbeck (OU) pro cess, such that 
the pattern evolves smoothly in space and time ( Schm idt et al.1 
120091: [Federrath et alll2010l) . The driving, and the PHANTOM 
Smoothed Particle Hyd rodynamics (SPH) code em ployed, are 
desc ribed in detail in iPrice & Federrafhl d2010l) (PF10, see 
also iFederrath et al.l 120101) . Calculations were evolved for 
10 dynamical times [defined as td = L/(2A4c s )], using 
only results after ltd such that turbulence is fully established 
(e.g. IFederrath et al.1 120091) . The amplitude of the driving 
force was adjusted to give RMS Mach numbers in the range 
1 < M < 20 by varying the energy input per Fourier mode 
proportional to the Mach number squared, i.e., E st i r oc M 2 
while the correlation time for the OU process was set to td 
(for the nominally input M). 

Most importantly, unless otherwise specified we have 
driven the turbulence using purely solenoidal Fourier modes. 
Thus, according to the heuristic theory of IFederrath et al.l 
(2010) we should expect a relationship of the standard form 
with b » 1/3. 

2.3. Measuring the density variance 

We consider a range of methods for measuring the density 
variance from the simulations. 



i) Measure the linear variance, a 2 ^- 



directly - with no 
assumptions about log-normality or otherwise - and fit 
the measured relation as a function of A4. 



a-M. relation in supersonic turbulence 
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> 512 3 parts -> AMR grid (eff. 8192 3 ) 
. 256 3 parts -> AMR grid (eff. 4096 3 ) 
3 parts -> 512 3 grid 

3 parts -> AMR grid, PDF fitted around mean 
3 parts -> AMR grid, PDF fitted around mean 
3 parts -> 512 3 grid, PDF fitted around mean 



> 512 3 parts -> AMR grid (eff. 8192 3 ) 
. 25 6 3 parts -> AMR grid (eff. 4096 3 ) 
- 128 3 parts -> 5 12 3 grid 
. 512 3 parts -> 5 12 3 grid 
, 512 3 parts -> 256- 1 grid 
] 256-' parts -> 256- 1 grid 
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FIG. 1. — Measured relationship between the (volume-weighted) standard 
deviation of the logarithm of density a s as a function of RMS Mach num- 
ber from a series of solenoidally-driven supersonic turbulence calculations. 
The points show time averages, with error bars showing (temporal) la de- 
viations. The dashed lines show the standard relation (Eq.fSJ with b = 1/3 
an d b = 1/2 while the dotted line shows the best fitting relationship found 
by Lemaster & Stone (2008) (Eq. f3j- Differences between directly measur- 
ing <j s (open circles, filled circles, and plus signs) compared to fitting the 
PDF around the mean (*,x and squares) are not significant (i.e., smaller 
than the time-dependent fluctuations). Overall, the results are consistent with 
b = 1/3, as expected for solenoidally-driven turbulence from Feder rath et al.l 
1 2008, 2010), and indistinguishable from the |LS08| best fit. 

ii) Measure the logarithmic variance a 2 s directly and fit the 
measured relation. Infer cr p /p assuming a log-normal 
PDF via Eq. ©. 

iii) Measure s and fit the measured relation. Infer a s using 
Eq. (]5]l and in turn a p /p using Eq. ©. 

iv) Determine the value of <j s that gives the best fitting PDF 
in a restricted range around the mean. Infer <J p /p using 
Eq. ©. 

The objection to method i) is that it is sensitive to the tails 
of the density distribution, where time-dependent fluctua- 
tions and i ntermittency effects can cause deviations from log - 
normality dKritsuket al.ll2007t iFederrath et ail l2010i IPFIO) . 
On the other hand no assumptions are made regarding the 
PDF, whilst methods ii)-iv) assume a priori that the PDF is 
log-normal, though for methods ii) and iii) only for obtaining 
the linear variance. Method iv) is the usual approach used to 
fit b for a given Mach number, if one additionally assumes the 
relationship given by Eq. (fJJ — an assumption we do not need 
to make here since a range of Mach numbers are examined. 

While the results are discussed in more detail below, es- 
sentially we find that methods ii)-iv) all give similar results 
for cr s , independent of numerical resolution, but that di- 
rect measurements of a p /p [method i)] are highly resolution- 
dependent. 

Volume-weighted variances were computed from the (mass 
weighted) SPH data by interpolating the density field to a grid. 
We found that this procedure gave better results than the di- 
rect calculation from the particles we have previously advo- 
cated (PF10), particularly at high Mach numbers (A4 > 10) 
where assuming that the volume element m/pis constant over 
the smoothing radius is an increasingly poor approximation. 
However, capturing the full resolution in the density field was 
found to require an adaptive rather than fixed mesh. 



FIG. 2. — Relationship between the linear and logarithmic density variance 
as a function of both intrinsic SPH resolution (number of particles) and the 
grid size used to compute the variances (see legend). While the measurements 
of o-g are resolution independent, there is a strong dependence on both the 
SPH and grid resolution in the directly measured linear variance, o^ip- Using 

an AMR grid to compute volume-weighted variances captures the full density 
field resolution in the SPH simulations, but even in the highest resolution 
calculations (512 3 particles), ^ p jp is severely underestimated compared to 
the expected exponential relationship (Eq.[6] dashed line) for M > 5. 

All plots show volume-weighted quantities, time-averaged 
over 81 snapshots evenly spaced between t/td = 2 and 
t/td = 10, with error bars showing the (temporal) ler devi- 
ations from these values. 

3. DENSITY VARIANCE - MACH NUMBER RELATION IN 
SUPERSONIC, ISOTHERMAL TURBULENCE 

3.1. a s as a function of Ad 

The direct measurements of the standard deviation in the 
log density, cr s , are shown in Fig. Q] from the results of calcu- 
lations performed using 128 3 , 256 3 and 512 3 SPH particles, 
using methods ii) and iv) (see legend). Dashed lines show the 
standard relation (Eq. [2]) with 6 = 1/3 and b = 1/2, whi lst the 
dotted line shows the best fitting relationship found by lLS08l 
(Eq.[3]l. Both the 6=1/3 curve, expected fo r solenoi dally- 
driven turbulence (IFederrath et alJl2008l l20Toh and the lLS08l 
fit show reasonable fits to the data, and indeed cannot be dis- 
tinguished given the time variability present in the calcula- 
tions. However, adopting b = 0.5 is clearly not consistent 
with our solenoidally-driven res ults. T he results are also con- 
sistent with our earlier findings (PF10), that showed a conver- 
gence in both grid and SPH methods towards b « 0.35 — 0.4 
at Mach 10. 

The measured value of a s cannot be compared to observa- 
tions, since only the 3D variance in the li near density can be 
observationally inferred (e.g. using the BFP method). Thus 
it is necessary to either measure or infer the linear variance 
from simulations to make this comparison. 

3.2. Direct measurement of a p/p 

A direct measurement of the linear density variance is dif- 
ficult even with high resolution simulations, demonstrated in 
Fig. |2] that shows the directly measured cr^,- as a function of 

<7g. The dashed line shows the expected relationship for a 
log-normal distribution (Eq.[6]). The exponential relationship 
is resolved only at the lowest Mach numbers, A4 < 5, cor- 
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FIG . 3 . — Directly measured or inferred (see legend) standard deviation of the linear den sity *T p /p as a function of RMS Mach number from the solenoidally- 
driven supersonic turbulence calculations. For comparison, observational determinations by Padoan et al. 1 1997a) and Brunt (2010) in IC5 14 6 and Taurus (respec- 
tively) are shown, together with the expected 6=1/3 and 6 = 1 linear relationships for solenoidal and compressive forcing (respectively) {Federrath et aU2008l 
120101) . including the corresponding data points from Federrath et al. (2010) (1024 3 grid; cyan triangles). Direct measurements of & p /p are resolution-limited 
(see Fig. [2}, although the values inferred by assuming Eq. (6) are upper limits, whereas the observations are likely to be lower limits. The discrepancy between 
solenoidally-driven simulations and observations indicates that some amount of gravity and/or compressive driving is necessary to explain the observational 
results. 



responding to o 2 s < 1.4, while at higher Mach numbers a p /p 
is a strong function of resolution, most easily demonstrated 
by interpolating the highest resolution SPH calculations (512 3 
particles) to fixed grids of decreasing resolution (see legend). 
This dependence is the reason why we eventually interpolated 
the SPH data onto an adaptive mesh, refined such that Ax < h 
for all cells within the smoothing radius 2h of any given parti- 
cle, in order to obtain a result for <J p /p that captures the max- 
imum resolution available in the SPH simulations, though we 
remai n limited by the intrinsic resolution of the simulations. 
IPF 101 found that SPH simulations at Mach 10 resolved a max- 
imum density at 128 3 particles similar to that captured on a 
fixed grid at 512 3 grid cells. This is consistent with the results 
here, where it is necessary to refine the grid to an effective 
8192 3 for the Mach 20 calculations employing 512 3 SPH par- 
ticles. It is also evident that fully resolving the strong fluctua- 
tions in the linear density at high Mach number is intractable 
with current computational resources. Our findings also sug- 
gest that the linear density variance is likely to be severely 
underestimated by limited observational resolution, so the re- 
sults in Taurus and I C5146 are almost certainly lower limits 
(see also lBrunlll20Tol) . 

3.3. <7p/p as a function of M: Comparison to observations 

The direct - but resolution limited - measurements for 
crp/p, computed without assumptions from the AMR grid, are 
shown in Fig.[3](see legend) [i.e., method i), above]. A better 
approach is to use the fact that the measurements for a s are 
resolution independent (Fig.[T), meaning that we can use the 
assumption of log-normality to infer the fully resolved value 
for cTp/p, i.e., using Eq. © (method ii). The standard devia- 
tions computed in this way are also shown, and as expected 



are consistent with a linear Opi-p-M- relation w ith 6 = 1/3 for 
sol enoida l forcing ((Federrath et al. 2008, 2010), and also with 
the LS08 best fit, the latter translated in terms of a p /p using 
Eq.©. 

We are now in a position to compare with the observational 
results. It should first be noted that the assumption of the 
cr s -tjp/p relation for the simulations essentially gives us an 
upper limit on cr p /p, whereas (see above) the finite resolution 
of the observations almost certainly gives a lowe r limit on the 
variance. The results from lPadoan et al.| (1 199 7a) (in IC5146: 
b = 0.5 ± 0.05 at Mach 10) and lBruntl (120101) (in Taurus, 
-Mrms = 17.6 ± 1.8 and a p/ p = 8.49±^;ff) are also plot- 
ted, both of which are consistent with 6 = 0.5 but clearly in- 
consistent with the calculations of purely solenoidally-driven 
turbulence. 

The other extreme is given by the b = 1 line in Fi g. [3] 
corres pondi ng to purely compressive forcing ( Fe derrath et al.l 
120081 1201 Oh . Data points from two 1024 3 grid simula- 
tio ns of purely so l enoidal and purely compressive forcing 
by Federra th et al.l (120081 120101) are also shown (cyan trian- 
gles). Clearly, purely solenoidal and purely compressive forc- 
ing seem inconsistent with the observations, while a mixture 
and/or the addition of gravity can fit the observations, best fit 
by a linear relation with b « 0.5. The o- p -M. plane, however, 
needs to be populated with many more observational mea- 
sures to draw more definite conclusions about, e.g., regional 
and evolutionary variations. 

4. DENSITY VARIANCE - MACH NUMBER RELATION IN MHD 
TURBULENCE 

Since molecular clouds are observed to be magnetised, we 
have additionally computed a series of Magnetohydrodynam- 
ics (MHD) calculations, driven identically to their hydrody- 
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FIG. 4. — As in Fig.ff]but for a series of 256 3 grid-based MHD calculations 
with field strength characterised by the ratio of gas-to-magnetic pressure [3 
(see legend). There is a general decrease in the measured variance in the 
MHD simulations at high Mach number, though no cle ar tren d with magnetic 
field strength. The best-fitting relationship found by LS08 for strong field 
MHD calculations (dotted line) is consistent with our ft = 1 results in a 
similar parameter range (M. < 6), but too steep at higher M. The < 0.05 
points refer to calculations employing (3 = 0.05, 0.01 and 0.02 at Mach 4, 
10 and 20 respectively. 



namic counterparts but with an initially uniform magnetic 
field threading the box. These have been com puted on a fixed 
grid u s ing th e fl ash code, as descr ibed in iFederrath et al.l 
d2010l) . lPF10l and lBrunt et ail d2010allbh (note that the driving 
routines are implemented identically in both the SPH and grid 
code). The magnetic field strength in the MHD calculations 
is characterised by the ratio of gas-to-magnetic pressure f3 = 



P/Pr 



in the initial conditions, where P„ 



/Mo- 



As previously we have also examined the effect of resolution, 
with a finding similar to the hydrodynamic case - namely that 
direct measurements of a p /p are strongly resolution affected 
(underestimated) - even at modest Mach numbers, similar to 
the results shown in Fig. [2] but that measurements of a s are 
resolution independent (at > 256 3 grid cells). 

Fig.|4]shows the results, similar to Fig.[T]for the hydrody- 
namic case and over-plotted with the hydrodynamic 6=1/3 



and 6=1/2 relationships (dashed l ines) a s well as the best- 
fitting MHD relationship found by ILS08L The most obvi- 
ous difference with MHD is that the density variances are 
significantly lower than their hydrodynamic counterparts at 
high (sonic) Mach number M. > 10 - though conversely 
marginally higher at lower Mach numbers. On the whole 
increasing the field strength seems to decrease the mean a s 
slightly. Whilst a complete MHD study is beyond the scope 
of this paper, the results clearly illustrate a shallower relation- 
ship than the hydrodynamic 6=1/3 curve at high Mach num- 
ber, with u a only weakly dependent on M in this regime (for 
(3 < 1). Thus, if anything, adding magnetic fields decreases 
the variance in the density field, worsening the discrepancy 
with observations. 

5. CONCLUSIONS 

We have measured the relationship between the density 
variance and the mean square Mach number from a series 
of simulations of supersonic, isothermal, solenoidally-driven 
turbulence over a wide range of Mach numbers (1 < Ai < 
20). We find that the standard relationship given by Eq. (O 
with 6 = 1/3 provides a good fit to the data over this 
range, consistent with the heuristic theory of Federra th et alj 
(2010) for solenoidal driving and with similar measurements 
by iLemaster & Stonel (120081) at lower Mach numbers. While 
it is difficult to measure the variance in linear density directly 
from simulations with finite resolution, the inferred relation- 
ship (Eq. Q] with 6 = 1/3) appears inconsistent with ob- 
servational determinations (6 w 0.5) in Taurus and IC5146, 
suggesting that additional physics such as gravity is impor- 
tant in these clouds and/or that some form of compressive 
driving is relevan t . Thi s is consistent with the fin dings of 
iKainulainen et al.l (120091) and lLombardi etal.1 (120101) for Tau- 
rus, where self-gravity is invoked to explain the deviation 
from log-normality in the high density tail of the (column) 
density PDF. Magnetic fields do not help to explain the dis- 
crepancy. 
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